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Abstract 

We have measured the cross section o"(e^e~ t:^t:~^['^)) at DA<I>NE, the Frascati 
(;/)-factory, using events with initial state radiation photons emitted at smaU angle 
and inclusive of final state radiation. We present the analysis of a new data set 
corresponding to an integrated luminosity of 240 pb~^. We have achieved a reduced 
systematic uncertainty with respect to previously published KLOE results. From 
the cross section we obtain the pion form factor and the contribution to the muon 
magnetic anomaly from two-pion states in the mass range 0.592 < Mtttt < 0.975 
GeV. For the latter we find A'^'^a^=(387.2±0.5stat±2.4exp±2.3th) x 10-^°. 
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1 Introduction 



The muon magnetic anomaly, a^, lias been recently measured at Brookhaven with 
an accuracy of 0.54 ppm [1]. The value of in the standard model, is found to 
differ from the experimental value by 2.8 to 3.4 standard deviations p|3] . The main 
source of uncertainty in the estimate of is the hadronic contribution, which is not 
calculable in perturbative QCD. The hadronic contribution, at lowest order, A'^'^°a^, 
is obtained from a dispersive integral over the cross section for e^e~— >■ hadrons 
[His] . The e"'"e~-^ TT+TT" channel accounts for ~ 70% of ^'^°a^ and ~ 60% of its 
uncertainty. 

It should be noted that the physically measurable cross section for e^e~^ vr+vr" , as 
such, cannot be used in the dispersive integral for two reasons. The first, obviously, 
is that the measured cross section is affected by initial state radiation (ISR) which 
must not be included in the contribution to the muon anomaly. Even the energy 
at the 7r7r7 vertex is different from the nominal e^e~ collision energy. The second 
reason is more a question of tradition and book keeping. The photon at the 7r"'~7r~7 
vertex, at lowest order, is a bare photon, i.e. without vacuum polarization. The 
measured cross section must be therefore corrected for both effects, as we discuss 
later. Final state radiation (FSR) from the pions must instead be included. The 
measured quantities therefore require corrections for the photon vacuum polariza- 
tion, for ISR, and to ensure that pion FSR is included, since some of the events 
with FSR might have been rejected. In our measurement there are some additional 
corrections, mostly due to ambiguities between ISR and FSR because we measure 
the dipion mass and not the e^e~ collision energy. 

In 2005, we published [6] a measurement of the dipion contribution A'^'^a^, using 
the method described in Sec. |21 using data collected in 2001 for /Ldt=140 pb~^, 
with a fractional systematic error of 1.3%. We discuss in the following a new and 
more accurate measurement of the same quantity (additional information can be 
found in [7]). 



2 Measurement of e'^e tt+tt cross section at DA$NE 

The KLOE detector operates at DA$NE, the Frascati 0-factory, a "small angle" 
e^e~ collider running mainly at a center of mass energy equal to the meson 
mass, iy~1020 MeV. At DA$NE, we measure the differential cross section for 
e^e~ — > 7r+7r~7 as a function of the tt+tt" invariant mass, M,r7r, for ISR events, and 
obtain the dipion cross section cTtttt = cr{e^e~ tt+tt") from [S]: 



dcrfee — > TXTx-^) 
s - 



dML 



= aUMl)H{Ml,s). (1) 

ISR 



Eq. [T] defines H, the "radiator function" . H can be obtained from QED calculations 
and depends on the e~^e~ center-of-mass energy squared s. In Eq. [T]we neglect FSR, 
which however is included in our analysis. The cross section we obtain is inclusive 
of all radiation in the final state. 
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3 Selection of e~^e vr+vr 7 events and background rejection 

KLOE collected a nominal integrated luminosity of about 2.5 fb~^ up to the year 
2006. The results presented here use data with /Ldt=240 pb~^ taken in 2002. The 
statistical fractional error on the 7t~^tt~ contribution to the muon anomaly, A'^'^a^, 
for dipion masses between 0.592 and 0.975 GeV is smaller than 0.2%. 

3.1 The KLOE detector 

The KLOE detector consists of a cylindrical drift chamber (DC) [9] and an electro- 
magnetic calorimeter (EMC) [10]. The DC has a momentum resolution of (Jp^/pi_ ~ 
0.4% for tracks with polar angle 9 > 45° . Track points are measured in the DC 
with a resolution in r-0 of ~ 0.15 mm and ~ 2 mm in z. The EMC has an en- 
ergy resolution of aE/E ~ 5.7%/ ^E (GeV) and an excellent time resolution of 

at ~ 54 ps / ^JE~{GeV) (B 100 ps. Calorimeter clusters are reconstructed grouping 
together energy deposits close in space and time. A superconducting coil provides 
an axial magnetic field of 0.52 T along the bisector of the colliding beam directions. 
The bisector is taken as the z axis of our coordinate system. The x axis is horizon- 
tal, pointing to the center of the collider rings and the y axis is vertical, directed 
upwards. A cross section of the detector in the y, z plane is shown in Fig. [TJ Events 
used in this analysis were triggered |TT] by two energy deposits larger than 50 MeV 
in two sectors of the barrel calorimeter. 

The 7r"'"7r~7 cross section diverges at small photon angle as (1 — cos^6'^)~^ oc 1/6'^ 
making FSR 7i~^7i~'j processes and (j) decays relatively unimportant. For this reason 
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we measure da/dM^^ at small photon angle. 

Figure [1] shows the fiducial volumes we use in the following for pions and photons. 
Note that the photon is not detected. Below we list the requirements for event 
selection. 

(1) Two tracks of opposite sign curvature, crossing a cylinder of radius 8 cm and 
length 15 cm centered at the interaction point, must satisfy 50° < 6 < 130° and 
Pj_ > 160 MeV or \pz\ > 90 MeV, to ensure good reconstruction and efficiency. 

(2) The (unobserved) photon direction, reconstructed from the two tracks above 
as = — (p+ + P-) must satisfy | cos^^| > cos(7r/12) (15°). In the following 
(^nn = n — 9^ is the polar angle of the dipion system. 

(3) The main background processes are e~^e~ —>■ e"'"e~7, fi'^fi~'j and (p ir^ir^Tr^ 
decays. Signal events are distinguished from e~^e~ —>■ e~^e~'j events by particle 
identification (PID), using a pseudo-likelihood estimator [12] for each track, 
L±, based on time of fiight, energy and shape of the cluster associated to 
the track. Electrons deposit most of their energy near the entrance of the 
calorimeter while muons and pions lose energy almost uniformly along the 
depth of the calorimeter. Events with both tracks having L± < 0, as in the 
lower left rectangle of Fig. [21 left, are identified as electrons and rejected. The 
efficiency for this selection is larger than 99.95%, evaluated on 7r+7r"7 samples, 
obtained from both data and Monte Carlo. The probability for the e+e~7 
events to be misidentified by the estimator as 7r"'"7r~7 is 3%. 

(4) The event must satisfy a cut on the track mass variable, Mtrk- Assuming the 
presence of an unobserved photon and that the tracks belong to particles of 
the same mass, Mtrk is computed from energy and momentum conservation: 

(v^ - V^|p+|2 + - \/|P-P + Kk)' - (P+ + P-)' = 

where p± is the measured momentum of the positive (negative) particle, and 
only one of the four solutions is physical. A cut is applied in the Mtrk^M^^ 
plane, as shown in Fig. [21 The requirement Mtrk > 130 MeV rejects fi'^fi~'-f, 
and further suppresses the fraction of e"'"e~7 surviving the pseudo-likelihood 
selection. While the M^^ dependent curve rejects tt^tt'tt^. 

About 3 X 10^ events pass these criteria. 

Residual e"'"e~7 , fi~^fi~'-f and tt+tt^tt'^ background contamination levels are eval- 
uated by fitting the Mtrk spectrum of the accepted events with a superposition 
of Monte Carlo (MC) distributions describing signal and backgrounds. The nor- 
malized contributions from signal and backgrounds are free parameters of the fit, 
performed for 30 intervals in M^^ of 0.02 GeV^ width. We also estimate the con- 
tribution from e^e~^ 7r'^7r~e"'"e~events, using the MC generator Ekhara [I3]. The 
errors on the background estimates refiect the uncertainty on the production mech- 
anism and the errors on the normalization coefficients from the fit. The results are 
given in Table [H Other possible backgrounds from [13] (/q, a)'y — > 7r"'"7r~7, 
e'^e~ —>■ p^TT^ 7r+7r~7, iv^v')! the radiative return to the u, e'^e~ — > 
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|^7iSR TT+vr 7r°7, are negligible after the acceptance cuts. Background from the 
process e~^e~ — > e'^e~e~^e~ is estimated to be well below 0.1%. 

3.2 Improvements with respect to the published analysis 

With respect to our published result [6], the present analysis profits from lower 
machine background and more stable DA$NE operation in 2002. Improved data 
filters were also developed. 

(1) A new trigger (L3) implemented at the end of 2001 reduced the loss of events 
rejected as cosmic rays events from 30% to 0.2%. The event loss is determined 
from a downscaled control sample taken without the enforcement of the cosmic- 
ray veto. 

(2) An offline background filter efficiency of 95% resulted in a large systematic 
uncertainty in our previous result. A new filter with 98.5% efficiency and neg- 
ligible systematic uncertainty has been implemented. A downscaled sample is 
retained to evaluate the filter efficiency. 

(3) We no longer require that the two pion tracks form a vertex. The uncertainty 
in the corresponding efficiency in the 2005 analysis is now removed. 

Finally, the Bhabha cross section has been reevaluated with better accuracy|15] 
than at the time of our previous measurement, as discussed below. 

3.3 Efficiencies, acceptance and systematic errors 

An improved simulation of the detector response allows determination of the ef- 
ficiency for the event-selection criteria described above. The generator Phokhara, 
including next-to-leading-order ISR [16] and leading-order FSR corrections, as well 
as simultaneous emission of one ISR and one FSR photon [17] has been inserted in 
the standard KLOE MC Geanfi [18]. 

Corrections are needed for the trigger and tracking efficiency. We compare MC 
efficiencies with efficiencies obtained from a data control sample, and correct the 
MC where discrepancies are present. 




Fig. 2. Left: PID estimator of the positive track vs. PID estimator of the negative track. 
Right: signal and background distributions in the M^j-i^-M^^ plane. 
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Trigger. The efficiency is obtained from a subsample of 7r+7r~7 events in which a 
single pion satisfies the trigger requirements. Then, the trigger response for the other 
pion is parametrized as a function of its momentum and direction. The efficiency 
as a function of M^^^^ is obtained using the MC event distribution and ranges from 
96% to 99%. The result is checked with a subset of the same sample, selected with 
an essentially independent drift chamber trigger, and evaluating the efficiency satis- 
fying the calorimeter trigger directly as a function of Mtttt. The constant fractional 
difference of 0.1% is taken as the systematic uncertainty. 

Tracking. The tracking efficiency for single pions is evaluated for each charge as a 
function of momentum and polar angle, using 0— > tt^tt~ti^ and 7r"'"7r~7 events, both 
identified on the basis of calorimeter information and the observation of a pion 
track. The two control samples are complementary: 0— tt'^tt^tt^ decays are more 
abundant, but do not cover the whole momentum range of interest. The efficiency is 
~ 98% and constant in M^^. The fractional difference between the results obtained 
with the two samples is 0.3%, which is taken as systematic error. 

Pion ID. Each track is extrapolated to the calorimeter and at least one cluster 
is searched for within a sphere of radius |rcnt — i"ciu| < 90 cm. rent and v^xa are the 
coordinates of the track entry point and of the cluster centroid. The value of 90 
cm is chosen to include pion fragments and minimizes the systematic error on the 
trigger efficiency correction. The efficiency is evaluated on 7r"''7r~7 events where a 
single track has an associated cluster identified as pion (L > 0), and parametrizing 
the probability of the other track to find an associated cluster with L > 0. This 
probability is larger than 99.9%. No difference on the pion cluster efficiency is found 
by varying the association radius; the systematic error is taken as negligible. 

Efficiencies for Mtrk cuts and acceptance are evaluated from MC, corrected to re- 
produce data distributions. The efficiency of the selection is shown in Fig. [31 

E — 1 

0.6 r 
0.5 r 
0.4 r 
0.3 r 
0.2 r 
0.1 - 

t 

0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

M^^^ [GeV^] 

Fig. 3. Global efficiency of 7r7r7 event selection. 
Systematic uncertainties are obtained as follows. 

• The systematic uncertainty due to the Mtrk cut is obtained by moving the cuts in 
the Mtrk^Af^^ plane within reasonable limit of the mass resolution and evaluating 
the difference in the 7r"'"7r~7 spectrum. We find a fractional difference of 0.2% 
(constant in M^^^^) which we take as systematic error. 

• Systematic effects due to polar angle requirements for the pions, 50° < 9 < 
130°, and of dipion, | cos^^^rl > cos(15°), are estimated by varying the angular 
acceptance by 1° around the nominal value. The systematic error decreases from 
0.6% at Ml^ = 0.35 GeV^ to 0.1% at M^^ = 0.64 GeV^ and is negligible for larger 
values (see Table [2]). 
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3-4 Luminosity 



The absolute normalization of the data sample is obtained from the yield of large 
angle (55° <9 < 125°) Bhabha-scattering events. The integrated luminosity, /Ldt, 
is obtained by dividing the observed number of these events by the effective cross 
section evaluated by the MC generator Babayaga [TH], inserted in the KLOE MC 
Geanfi [20] • The Babayaga generator includes QED radiative corrections via the 
parton shower algorithm. An updated version of the generator, Babayaga@NLO 
[TSj . gives a Bhabha cross section of 456.2 nb, 0.7% lower than the value from the 
previous version, while the theoretical uncertainty is reduced from 0.5% to 0.1%. The 
experimental uncertainty on the luminosity is 0.3%, dominated by the systematics 
on the angular acceptance. 

4 Evaluation of the pion form factor 

The differential 7r+7r~7 cross section is obtained from the observed number of events, 
A'obs, after subtracting the residual background, A^bkg, correcting for the selection 
efficiency, esei(M^^), and dividing by the luminosity L: 



The background mentioned above varies smoothly from 1% around the p peak to 
~7%, for low and high M^^ mass. This background is dominated by misidentified 
muon pairs at high mass, and 7r"'"7r~7r° events at low mass. Our mass resolution 
{5Ml^ ~ 2 X 10-3 GeV^ or (5M^^~1.3 MeV at the p peak) allows us to use bins of 
AM^^ = 0.01 GeV^ width. In order to correct for resolution effects, the differential 
cross section is unfolded using the Bayesian method described in [21]. The point-by- 
point uncertainty introduced by the unfolding procedure - appreciable only in the 
p-uo region, M^^ ~ 0.6 GeV^ - is given in Table [3l The unfolding does not introduce 
any additional systematic error on A'^'^a^. 

The quantity M^.^ is computed from measured momenta of the pions and is shifted 
by radiative effects from the mass value at the tt^'k~'^ vertex, M°^. The cross 
section (Jt^.,^{M^^) is determined accounting for this shift and dividing the ti^-k^'-) 
cross section by the radiator function H (obtained from Phokhara by setting pion 
form factor = 1) as in Eq. [TJ Fluctuations in the e^e~ CM energy during data 
taking introduce an additional systematic uncertainty in the evaluation of H (see 
Table H. 

The pion form factor is obtained from ct^tt after subtraction of FSR, assuming point- 
like pions {r]FSR term [22|23] ): 





1 



(2) 



i^.|'(l+r/FSR) 




3 Ml 




(3) 



a- 



TTTT •> 



'vr 



\ M2 

\ ^-^ TTTT 
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where a is the fine structure constant [a = e^/Air), and [3.^ is the pion velocity in 
the vrvr rest frame. 

Our results are summarized in Table [71 which gives the integrals over bins of 0.01 
GeV^ identified by their lower edge of: 

• the observed differential cross section dae+e-^nwy/dM^^ as a function of the in- 
variant mass of the dipion system in the angular region Ot^tt (tt — 6-,^^) < 15°; 

• the bare cross section cr°^, inclusive of FSR, and with vacuum polarization effects 
removed [21]: a^^ = a^^ [«(0)/a(M,,)]'; 

• the pion form factor without FSR and with vacuum polarization effects included. 

The errors in Table [7] are statistical only. The systematic errors are given in Ta- 
bles [H [31 and [H The statistical errors are weakly correlated as a result of the 
resolution unfolding. The correct covariance matrix is used for the calculation of 
the statistical error on A'^'^a^. The systematic errors cannot be considered as fluc- 
tuations of random variables within each of the categories considered. We combine 
all contributions for the same M^^ value in quadrature. We add linearly then the 
errors for each bin to obtain the total systematic error on A'^'^a^. 

Figure [H left, shows the differential cross section for e'^e~ ^ 7r~^ it ~'~f with | cos6'^| > 
cos(15°) after applying the corrections described above. Figure [H right, shows the 
cross section a^^, which is the input for the dispersive integral for A'^'^a^. 
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Fig. 4. Left: Differential cross section for e~^e 
Right: cross section for e"'"e~ — > 7r~^7r~. 



•vr+vr 7(7), with |cos0^| > cos(15°). 



5 Evaluation of A'^'^a^ 

The dispersive integral for A'^'^a^ is computed as the sum of the values for listed 
in Table [7] times the kernel K{s): 

^^^""^^i^ / , (4) 
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Table 1 

Systematic errors in (in percent) due to background subtraction, in 0.01 GeV^ bin intervals 
of M^^. The lower bin's edge is given by the sum of the values in the first row and first 
column. 



M'^^ range ( GeV^) 


Systematic error (%) 
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Table 2 

Fractional systematic errors on the acceptance. 
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Table 3 

Systematic error (in percent) on dag+g-^^^^/dM^ due to the correction for detector 
resolution in 0.01 GeV^ intervals. The indicated values for represent the lower bin's 
edge. Outside this interval, the p-uj region, the effect is negligible. 



(Ml)' (GeV2) 


0.00 


0.01 


0.02 


0.03 


0.04 


0.05 


0.06 


0.07 


0.08 


0.09 


0.3 












0.3 


0.4 


0.2 


0.1 


0.2 


0.4 


0.1 


0.2 


0.1 


0.2 


0.1 


0.2 


0.2 


0.3 


0.1 


0.2 


0.5 


0.3 


0.3 


0.1 


0.2 


0.3 


0.3 


0.3 


0.2 


0.2 


0.3 


0.6 


0.1 


0.3 


0.3 


0.2 


0.2 


0.3 


0.3 


0.3 


0.2 


0.3 


0.7 


0.3 


0.3 


0.4 


0.3 


0.3 


0.3 


0.4 


0.3 


0.3 


0.3 


0.8 


0.4 


0.4 


0.3 


0.5 


0.4 


0.5 


0.4 


0.4 


0.5 


0.5 


0.9 


0.5 


0.6 


0.6 


0.6 


0.7 













Table 4 

Systematic errors (in percent) on the radiator function due to the spread of ^/s in the 
2002 data taking period, given in 0.01 GeV^ bin intervals of M^. The lower bin's edge is 
given by the sum of the values in the first row and first column. 



where the kernel, see the second paper of ref. [S], is given by 

K{s) = (l - y ) + (1 + XY{1 + X"2) (log(l + X) - X + y ) + ^^-^X^ log X 
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Reconstruction Filter 
Background subtraction 
Trackmass 
Pion cluster ID 
Tracking efficiency 
Trigger efficiency 
Acceptance 
Unfolding 
L3 filter 

^/s dependence of H 
Luminosity 


negligible 


Tab.m 


0.3% 


0.2% 
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0.3% 


0.1% 


Tab. [2] 


0.2% 


Tab. El 
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0.1% 




Tab.H 


0.2% 


0.3% 


Experimental systematics 




0.6% 


FSR resummation 
Radiator function H 
Vacuum Polarization 




0.3% 




0.5% 




0.1% 




0.1% 


Theory systematics 




0.6% 



Table 5 

Systematic errors on 0"7r7r7, cr^J^^, F^^ and IS.'^'^a^. 
with 

_ 1 - y^l - 4m2/g 
l + 7l-4m2/s 

Eq. 1 gives A^^a^ = (387.2 ± O.S^tat ± 2.4exp ± 2.3th) x 10"^° in the interval 0.35< 
M^^ < and 0.95 GeV^. Contributions to the systematic errors on A'^'^a^ are given 
in the last column of Table [51 

6 Comparison between 2008 and 2005 analyses 

In order to compare consistently the tt^tt^'j differential cross section from this 
analysis to that from our previous analysis, two corrections have been applied to 
the previous results: 

• a —0.7% overall shift, due to the new evaluation of the Bhabha cross section, 
obtained from the updated version of the Babayaga generator (see Sec. 13. 4p : 

• an energy- dependent effect due to a double counting of the calorimeter cluster 
efficiency in the evaluation of the trigger correction, which overestimates the cross 
section mainly at low-mass values by a few percent. 

As a result of these updates, the value of A'^'^a^ from our previous analysis changes 
to (384.4 ± O.Sstat ± 4.6sys) x 10~^°. The fractional difference between the spectra 
for the present analysis and that previously published (with updates), is shown in 
Fig.[5l While the agreement below the p peak is good, above 0.7 GeV^ there is some 
difference between the two spectra. The value obtained for the integral is consistent 
between the two data sets (as shown in Table [6]). Because of the improvements to 
the analysis described in Sec. 13.21 and the high quality of the 2002 data, we consider 
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d(T(7r7r7, '08)/da{Tnrj, '05) — 1 




■ M^l_ (Gef) 

0.4 0.5 0.6 0.7 0.8 0.9 
Fig. 5. Comparison of the present result with the pubhshed data, updated for the effects 
described in the text. The band is just the fractional systematic error of the ratio. 

the present result to supersede that previously published. 



A'^'^a^.xlO^o 0.35 < M^^ < 0.95 GeV^ 


published 05 


388.7 ± 0.8stat ± 4.9sys 


updated 05 


384.4 ± 0.8stat ± 4.6sys 


new data 08 


387.2 lb O.Sstat i 3.3sys 



Table 6 

Comparison among A'^'^a^ values from KLOE analyses. 



7 Comparison with CMD-2 and SND results 

We may compare the present result on l-FVP with the results from the energy scan 
experiments at Novosibirsk CMD-2 [25] and SND|26]. For a given energy scan ex- 
periment, whenever there are several data points falling in one 0.01 GeV^ bin, we 
average the values. The result can be seen in Fig. [6l left. Figure [6], right, shows 
the fractional difference between the data points from the energy scan experiments 
(CMD-2 and SND) and the KLOE data. There is reasonable agreement between 



i CMD-2 , ' 

- I • SND M^„(GeV^) 

3 0.4 0.5 0.6 0.7 0.8 0.9 1 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Fig. 6. Left. from CMD-2 [25], SND [26] and KLOE with statistical errors. Right. 

Fractional difference between CMD-2 (C) or SND (S) and the KLOE (K) results. The 
dark (light) band is the KLOE statistical (statistical © systematic) error. Also for CMD-2 
and SND statistical © systematic errors are shown. 

the experiments, as also indicated by the computed values of A'^^'a^ given below in 
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the range of overlap 0.630< Afn-Tr < 0.958 GeV, combining statistical and systematic 
errors in quadrature: 

SND, 2006 [26] A^^a^, = (361.0 ± 5.1) x 10"^° 
CMD-2, 2007 A^^a^ = (361.5 ± 3.4) x 10"^° 
this work A^^a^ = (356.7 ± 3.1) x 10~^°. 

A fit for the best value gives 359.2±2.1 with x^/dof=1.24/2, corresponding to a 
confidence level of 54%. 

8 Conclusions 

We have measured the dipion contribution to the muon anomaly, A'^'^a^, in the 
interval 0.592 < Mj^^ < 0.975 GeV, with negligible statistical error and a 0.6% 
experimental systematic uncertainty. Radiative corrections calculations increase the 
systematic uncertainty to 0.9%. Combining all errors we find: 

A^X(0-592<M^^< 0.975 GeV) = (387.2 ± 3.3) x 10"^°. 

This result is consistent with our previous value, with a total error smaller by 30%. 
Our new result confirms the current disagreement between the standard model 
prediction for and the measured value. 
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Table 7 

da-jYTT-y/dM^^, cross section and the pion form factor, in 0.01 GeV^ intervals. The value 
given in the M^^ column indicates the bin center. 
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